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Abstract 

The aim of this study is to address the effects of wind-induced drift on a floating sea objects using high- 
resolution ocean forecast data and atmospheric data. Two applications of stochastic Leeway model for 
prediction of trajectories drift in the Mediterranean sea are presented: long-term simulation of sea drifters 
in the western Adriatic sea ( 21 . 06.2009 — 23 . 06 . 2009 ) and numerical reconstruction of the Elba accident 
( 21 . 06.2009 — 23 . 06 . 2009 ). Long-term simulations in the western Adriatic sea are performed using wind 
data from the European Center for Medium-Range Weather Forecast (ECMWF) and currents from the 
Adriatic Forecasting System (AFS). An algorithm of spatial clustering is proposed to identify the most 
probable search areas with a high density of drifters. The results are compared for different simulation 
scenarios using different categories of drifters and forcing fields. The reconstruction of sea object drift near 
to the Elba Island is performed using surface currents from the Mediterranean Forecasting System (MFS) 
and atmospheric forcing fields from the ECMWF. The results showed that draft-limited to an upper surface 
drifters more closely reproduced target trajectory during the accident. 
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1. Introduction 

The task of developing fast and reliable 
computer-assisted methods for search and rescue 
(SAR) operations becomes increasingly important 
due to increased sea traffic and higher risk of emer¬ 
gencies at sea pQ. 

In previous decades the existing SAR systems 
relied mostly on climatological means as well as 
were lacking high resolution environmental data 
mi- Advances in a high resolution operational 
forecast in the Mediterranean sea mmmm pro¬ 
vided a framework for a more accurate prediction 
of drifters forecast in this area. Also the inclusion 
of assimilation schemes in ocean forecast models [7] 
helped to improve a prediction of time-dependent 
currents and wind with respect to earlier models. 

A large databank of windage characteristics of 
sea objects was collected over the years that 


uncovered a linear relationship between drift veloc¬ 
ity and applied surface wind for different categories 
of sea drifters. The field observations were used in 
the context of trajectory prediction using the ap¬ 
proach of stochastic Leeway drift 10 ITT]. The Lee¬ 
way model m produced final positions of drifters 
based on their windage characteristics [10. Addi¬ 
tionally the modeling approach enabled prediction 
of a search area in terms of a probability density. 
However, the model does not classify the drifters 
according to their temporal and spatial behavior. 
While backtracking simulations could provide a full 
classification of spatially separated trajectories time 
limitation could be an issue in the SAR operations 
and a faster and simpler computational alternative 
should be provided. 

Here we proposed a computational approach for 
characterization of drifters based on a density clus¬ 
tering algorithm m • It was designed to improve 
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qualitative information about the stochastic out¬ 
come of the Leeway simulations and to increase the 
prediction probability. The spatial clustering local¬ 
izes the areas of higher probability of drifters and 
finds possible outliers. Also the procedure sorts out 
ensemble of drifters according to their temporal be¬ 
havior. 

We used the modeling framework based on the 
explicit parameterization of an object drift by wind 
P3B nn. Two applications of the stochastic Lee¬ 
way model mm for tracking of small size sur¬ 
face drifters were studied: (1) long-term simu¬ 
lation of sea drifters in the western Adriatic sea 
(21.06.2009 — 23.06.2009) and (2) numerical re¬ 
construction of the Elba accident (21.06.2009 — 
23.06.2009). In the first study experiments were 
performed with drifters released in the vicinity of 
the known hyperbolic structure of the flow near to 
the Gargano peninsula in the western Adriatic sea. 
The spatial sorting of drifters was studied by the 
determination of distinct clusters. Different cate¬ 
gories of drifters and geophysical forcing were tested 
as well. In the second study qualitative comparison 
between the behavior of reconstructed trajectories 
of drifters and known information about the Elba 
accident was performed using ensembles of stochas¬ 
tic simulations. In both applications the clustering 
method provided improved probability estimate of 
final search areas via the sorting of drifters accord¬ 
ing to their spatial distribution. 

2. Model 

2.1. Model description 

The Leeway drift is defined as a drift of floating 
object with respect to ambient current under the 
influence of wind and waves EDnsiini. According 
to the definition the leeway velocity is determined 
from a difference between drift velocity Vd r and 
Eulerian current velocity Ve of a sea object: 

L = V dr - V E , (1) 

where L is the Leeway velocity. The integration tra¬ 
jectory is obtained by calculating the drifters path 
according to formula: 

dr = Vdr dt T de — (L T Ve)^£ T de, (2) 

where dt is the time step, de represents a diffusive 
component of velocity due to sub-scale variability 
of ocean currents that is below the given model res¬ 
olution, Ve is the Eulerian current taken from the 


ocean model. The Leeway drift L is parametrized 
by 10 m atmospheric wind W\o (see Appendix 
Eq. [6|. Since in field campaigns the contribution 
of the Stokes drift to the motion of a drifter was 
difficult to separate from the effect of wind 0. In 
the present modeling paradigm a sea dominated 
by weak waves is considered and the waves are as¬ 
sumed to be alighed with the wind. The influence 
of Stokes drift is implicitly included through Lee¬ 
way coefficients (see Appendix Eq. [6|. The model 
accounts for the action of wind on small drifters and 
neglect damping and excitation by waves mm- 

Since in general shape and size of a sea object 
is not regular it moves at a certain angle to the 
wind direction m Therefore, the Leeway veloc¬ 
ity is decomposed into the downwind (DWL) and 
crosswind (CWL) Leeway components which are 
projections of Leeway drift on wind direction Wi 0 
mm- Also due to uncertainty of initial orienta¬ 
tion of a drifter with respect to wind an initial en¬ 
semble is generated with a half of drifters oriented 
left and remaining half right with respect to wind 
direction m ■ Because of the existing variability of 
wind, currents and some perturbation of the drifter 
motion its orientation can abruptly change. The 
change of orientation is introduced via the alterna¬ 
tion of the sign of the CWL component. In this 
study probability of change of sign equal to 4% per 
integration step j8]. 

Since a finite size of geographical domain was 
considered the drifters could reach the bound¬ 
ary of the domain or approach the coastal shore¬ 
line. In the former case they were marked as 
off-grid and removed from further consideration. 
Stranded drifters were identified according to a 
high-resolution coastline contour [18]. 

Often during emergency situations at sea exact 
information about an accident is missing. The Lee¬ 
way model accounts for the uncertainties in the po¬ 
sition of an accident m by assigning initial posi¬ 
tions from a normal distribution centered around an 
expected location of an accident. Particularly, the 
drifters are released continuously in a spatial area 
defined by initial distance ro = 0.5cr around the 
last known position (LKP), where a r is the vari¬ 
ance of drifters initial positions. The LKP refers to 
the possible location of an accident and is provided 
in terms of longitude and latitude for every drifter. 
In the experiments studied here drifters were re¬ 
leased simultaneously. In general uncertainty in 
time could be included m- 

To include the sub-grid variabilities of the wind 
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and current velocities the Gaussian perturbations 
were added at every integration step t n and for the 
fcth drifter: 

curr / -l r\ _ curr , ccurr ( o\ 

u x,y \ L n^) — V x ,y ' $x,y > W 

w x ,y{t n ,k) = w x>v +£™™ d , ( 4 ) 

where and ^y d are th e uncertainties cho¬ 

sen from a normal distributions Af(0,cr curr ) and 
N(0iCJwind) correspondingly. The standard devi¬ 
ation (std) of wind cr w i n d current cr curr were esti¬ 
mated using averages over the corresponding geo¬ 
graphical domain. 

3. Methods 

3.1. Interpolation schemes 

The Eulerian velocity of drifter Ve was ob¬ 
tained using bilinear interpolation of the model 
zonal and meridional components at the position 
of a drifter. Also weighted linear time interpolation 
was used to interpolate current and wind velocities 
between two subsequent time stamps. In our sim¬ 
ulations the integration step St = 360s was taken. 
For the advection scheme the Runge Kutta second- 
order method was used with the Euler first-order 
trial step 20. E]. 

3.2. Clustering procedure 

In the implementation of the algorithm the fi¬ 
nal positions of drifters were used for evaluation of 
a spatial separation between pairs of drifters. The 
separation distance was estimated from the haver- 
sine formula: 

d = 2 i ?0 arcsin ^Jdj + cos li cos IjhCj > ( 5 ) 

where R 0 is the Earth radius, i,j = 1... Nd r 
are drifter indices, A*,</>* are longitudes and lati¬ 
tudes of the fcth drifter and f Xl ,x 2 = ski 2 Xl ^ X2 . 
The mean separation distance D mean was evaluated 
from Eq. [5] across all distinct pairs of drifters. Also 
an upper bound D max > D mean was assigned to 
exclude mean separation distance that exceeds the 
size of a given geographical domain. Indeed, in the 
case Dmean > D ma x no clusters were searched. 

For the fcth drifter a new cluster was formed by 
the inclusion of the fcth neighbor provided that the 
condition D ^ < D me an was satisfied. The proce¬ 
dure was finalized after sorting all possible pairs of 
drifters. 


Application of the algorithm produced a set of 
spatially separated clusters. We classified a proba¬ 
bility of containment ( POC ) of drifters per cluster 
as follows: POC = N^/N^r x 100 %, where Nf, is 
a number of drifters inside the fcth cluster and Nd r 
is an ensemble size. The mean cluster trajectory 
(center of mass of the cluster) is defined as the time 
sequence of mean longitudes and latitudes across all 
cluster members calculated at every time step. 

4. Oceanographic and atmospheric data 

Six hourly 10 m wind fields with 0.5° x 0.5° spa¬ 
tial resolution were retrieved from the ECMWF. 

For typical SAR objects of small size uni only 
top surface levels of ocean currents were used. In 
our study, we used oceanographic data from two 
operational models: the MFS with a horizontal res¬ 
olution of about 1/16° x 1/16° [4[ [ 6 ] and the AFS 
(horizontal resolution from 1 / 22 ° to 1/45°) [5 ; T9] . 
In the vicinity of shorelines a sea-over-land inter¬ 
polation procedure was implemented to provide a 
horizontal interpolation of the ocean currents to¬ 
wards coastal zone. 

5. Case study I: experiments in the Adriatic 

sea 

The stochastic simulations of an ensemble of 
Ndr = 0 ( 2000 ) drifters were performed using var¬ 
ious simulation scenarios (see the summary of ex¬ 
periments in Table [l]). All drifters were released 
simultaneously from the position 42°35' N, 16° E. 

The mean surface wind evaluated for a study 
period (Fig. [l]) strongly contributed to the north¬ 
eastern displacement of surface drifters in the area 
around the release position. Also the western Adri¬ 
atic boundary current (Figs. [ 2 ]) was influencing the 
transport of near coast drifters southwards during 
the given period. 

In the first experiment (Fig. [3]) a significant part 
of an ensemble of drifters was carried by the western 
coastal flow initially southwards and later towards 
the inner southern Adriatic circulation. The major¬ 
ity of drifters were localized in the southern Adri¬ 
atic (clusters 2 and 4), a small part of the ensemble 
(cluster 1 ) was stranded on the western coast and 
the remaining drifters due to low velocities were lo¬ 
calized inside the release zone (cluster 3). 

The second experiment (Fig. [4| was done for 
person in water (PIW) category (see Table [ 2 | with 
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Table 1: Summary of parameters for the simulation cases: experiment number, category of drifter ((-) sign is indicated when 
no wind is used), initial radius ro of release around LKP [km] and depth [m]. 


Exp. 

Category 

ro 

Depth 

1 


10 

0 

2 

PIW (mean) 

10 

0 

3 

- 

10 

10 


the atmospheric wind. During 21 days of simula¬ 
tions the wind and currents influenced the disper¬ 
sion of drifters across the whole Adriatic basin to¬ 
wards the eastern Adriatic coast: small size clusters 
(clusters 2,4 and 5) and isolated drifters were local¬ 
ized in the central Adriatic. As in the experiment 
without the wind (Fig. [3| a fraction of stranded 
drifters was found inside the release area (cluster 
1). No large size clusters were formed in this case. 


30.10.2008-20.11.2008 



14°E 15°E 16°E 17°E 18°E 19°E 
mean wind speed, m/s 

V'V'b fc <b <o A % ^ ^ 

I I 

I Reference Vector | 


Figure 1: (Color online) (a) Magnitude, direction of mean 
[m/s] of the ECMWF winds for 30.10.2008 - 20.11.2008. 


30.10.2008-20.11.2008 

44°N 

43°N 

42°N 

41 °N 



14°E 15°E 16°E 17°E 18°E 19°E 


mean current speed, m/s 




Figure 2: (Color online) Magnitude and direction of the 
mean AFS currents for 30.10.2008 — 20.11.2008: (a) mean 
surface currents [m/s]. 


In the third experiment (Fig. [5| a current at 
the 10-m depth level was used to account for pos¬ 
sible scenarios of motion of submerged drifters. 
A higher clustering tendency was found in this 
case: the drifters were densely positioned inside 
three large clusters. The largest group of drifters 
(cluster 1 with POC > 50%) was localized in¬ 
side the southern Adriatic circulation. Also the 
western Adriatic current contributed towards the 
transport of a group of drifters (cluster 3) to 
the southern end of the Salento peninsula. No 
stranded drifters were found apart from a set 
of drifters that did not leave the release area. 



Figure 3: (Color online) Leeway simulations using the AFS 
daily surface currents without wind: the LKP (big dot), final 
position of drifters (red dots), mean trajectories of drift for 
every cluster (black curve), cluster numbers and search area 
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(a) 


(black convex polygon) for every cluster, total area of search 
with all drifters inside (black convex polygon), drifters- 
outliers in open sea (red dots) and stranded drifters (blue 
triangles). Four clusters were found: (1) POC = 8.55%, (2) 
POC = 41.45%, (3) POC = 5.65% and (4) POC = 22.85%. 



Figure 4: (Color online) Leeway simulations with inclusion 
of 10 m wind. Values for PIW (mean values) are used to 
parametrize the Leeway drift. Five clusters were formed: 
(1) POC = 17.1%, (2) POC = 6.9%, (3) POC = 5.9%, (4) 
POC = 6.6% and (5) POC = 3.6%. 



Figure 5: (Color online) Leeway simulations for subsurface 
drift. The current at ~ 10 m depth of the MFS model was 
used. Three clusters were formed: (1) POC = 50%, (2) 
POC = 42.6% and (3) POC = 4%. 


6. Case study II: reconstruction of the Elba 

accident 

To set up initial conditions and define param¬ 
eterization of an object drift the documented data 
about the accident were used [22]: a person on an 
inflatable raft was lost in vicinity of coast of the 
Elba Island at 42°43 / 60 // N, 10 o 9'5" E on 21.06.2009 
at 1 : 30 UTC and after 34 hours he was found 
alive at the position 42°22'90 // N, 9°53 / 30 // E in the 
open sea, it is known that at the time of rescue 
(22.06.2009 at 11 : 30 UTC) the boat was partially 
deflated. 


21.06.2009 
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42°40'N 
42°30’N 
42°20'N 
42° 1 0'N 
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wind speed, m/s 


Figure 6: (Color online) daily wind direction and magni¬ 
tude from the ECMWF for days: (a) 21.06.2009 and (b) 
22.06.2009. 

42°50'N 
42°40'N 
42°30'N 
42°20'N 
42°1 O' N 

9°40'E 10°E 10°20'E 


42°50' N 
42°40'N 
42°30'N 
42°20'N 
42°10'N 

9°40'E 10°E 10°20'E 

current speed, m/s 


current speed, m/s 

(b) 22.06.2009 
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Figure 7: (Color online) daily surface currents direction and 
magnitude from the MFS for days: (a) 21.06.2009 and (b) 
22.06.2009. 

Since the exact object parameterization was 
not known the simulations of drifters were per¬ 
formed using the categories similar to a PIW on 
an inflatable boat (see Table [2|. The ensemble of 
Ndr 0(3000) members was released simultane¬ 
ously from the same LKP (42°43 / 60 // N, 10°9 / 5 // E). 
The start date and time 30.10.2008,10 : 30 UTC of 
the release of drifters was the same in all the experi¬ 
ments. Drifters were released in a circle with radius 
ro = 15 km around LKP. The minimum number of 
drifters per cluster is N m i n = 20. 

Both the south-western wind on 21.06.2009 
(Fig. |6] a) and the westward surface current 
(Fig. [7pa) contributed to initial southwestern dis¬ 
placement of drifters with respect to the release area 

To summarize the results of simulations: the 
trajectories of drifters with higher DWL and lower 
CWL (see Table [2| more closely reproduced target 
trajectory in the accident. Generally, these drifters 
are draft-limited to an upper sea. Hence their 
trajectories were more influenced by the wind and 
wind-induced forcings during experiments. In con¬ 
trast, the vessels with a deep draft were less affected 
by wind and drifted westwards from the Elba loca¬ 
tion. Due to details of ensemble initialization a sig¬ 
nificant fraction of drifters was found stranded on 
the neighboring island and a set of drifters with low 
initial velocities was stranded near to the LKP. In 
the experiments with draft limited drifters (Figs. [8] 
d and e) a larger separation between sets of coastal 
and open sea drifters was observed. 

7. Conclusions 

The stochastic simulations in the Adriatic sea 
identified most probable search areas with a high 
concentration of drifters: while a significant portion 
of ensemble members entered the central Adriatic 
pit another group of drifters moved with the west¬ 
ern Adriatic coastal current and was either found on 
the coast or entered the southern Adriatic circula¬ 
tion. We observed high sensitivity of trajectories to 
the details of an initial ensemble such as category of 
drifters as well as choice of geophysical forcing. In 
all the experiments the central and southern Adri¬ 
atic circulations were important contributions for 


Figs. [8] On the following day the decrease of the 
mean surface current (Fig. [7] b) and change of the 
wind direction (Fig. [6] b) favored a slight deflec¬ 
tion of drift trajectories to the southeast (Figs. [8]). 
In the first experiment (Fig. [8] a) the search area 
included the found position but the drifters were 
spread over large area and thus their average den¬ 
sity around the found position was low. In the next 
two experiments (Fig. [8] b and c) an overall south¬ 
western drift from the LKP was observed. Two 
clusters were found for the liferaft with a deep bal¬ 
last with the largest (POC > 91%) formed by open 
sea drifters. In the last two experiments (Figs. [8] d 
and e) the southwards drift was more pronounced 
and trajectories more closely represented the ob¬ 
served real trajectory during the accident. The es¬ 
timated search areas (Figs. [8] d and e) included the 
position of rescue with the largest clusters formed. 

the transport of drifters. These circulations me¬ 
diated the detachement of drifters from the coast 
and further displacement inside the main Adriatic 
gyres. The clustering procedure helped finding very 
localized search areas with a high density of drifters 
in the experiment with subsurface drifters. Our 
results confirmed previous studies [12 j that high¬ 
lighted an influence of underlying flow hyperbolic- 
ity on the transport of drifters near to the Gargano 
Cape. 

In the numerical reconstruction of the Elba acci¬ 
dent the estimation of the most probable drift path 
was provided. The search area identified by cluster¬ 
ing method matched the documented final position. 
The results showed that draft-limited drifters more 
closely reproduced observed final position during 
the accident. 

8. Appendix 

8.1. Definitions of Leeway coefficients 

Here the explicit form of empirical relation¬ 
ship between Leeway components and magnitude 
of wind is provided. The DWL component Ld and 
CWL component L c are linearly related to the wind 
magnitude Wio = according to derived 

relationships mm- 

Ld,c ad,cWio T (6) 

where a^, bd and a c , b c are the linear regression co¬ 
efficients for the DWL and CWL correspondingly. 
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(a) 21.06.2009-23.06.2009 (b) 21.06.2009-23.06.2009 






(e) 21.06.2009-23.06.2009 



Figure 8: (Color online) Leeway simulations for the Elba accident: LKP (42°43 / 60 // N, 10°9 / 50 // E), position where the person 
was found (42°22 / 90 // N, 9°53 / 30 // E ), final positions (red dots) in open sea and (blue triangles) stranded on coast, area of 
drifters release (red convex polygon), final area of search (black convex polygon), (a) Results for the a life raft without ballast. 
Clustering algorithm distinguishes a single cluster with all drifters included, (b) Results for a life raft with a deep ballast. Two 
clusters are found: (1) POC = 7.62% and (2) POC = 92.1%. Cluster 1 consists of stranded drifters, (c) Results for a surfboard 
with person. Drifters form a single cluster, (d) Results of simulation for a sports boat. Three clusters are distinguished: (1) 
POC = 19.6% (stranded drifters), (2) POC = 7.2% (stranded drifters) and (3) POC = 72.8% (open sea drifters), (e) Results 
for a sport fisher. Three clusters are found: (1) POC = 73.1% (open sea drifters), (2) POC = 19.8% (stranded drifters) and 
(3) POC = 6.6% (drifters near to the LKP). 
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Table 2: Leeway coefficients for the categories used in the simulations. The regression slope values a^c [%], the offset 
[cm/s] and the standard deviations cr^ ;C [cm/s] are shown for the DWL and the left/right CWL components. Values are 
provided from [2]. 


Object name 

&d 

DWL 

Pd, 

°d 

right 

Oi c 

CWL 

Pc, 


left 

a c 

CWL 

Pc 


Life raft, no ballast 

3.7 

0 

12 

1.98 

0 

9.4 

-1.98 

0 

9.4 

Life raft, deep ballast 

3.52 

-2.5 

6.1 

0.62 

-3 

3.5 

-0.45 

-0.2 

3.6 

Surf board with person 

1.93 

0 

8.3 

0.51 

0 

6.7 

-0.51 

0 

6.7 

Sport-boat 

6.54 

0 

3.0 

2.19 

0 

2.8 

-2.19 

0 

2.8 

Sport-fisher 

5.55 

0 

3.3 

2.27 

0 

3 

-2.27 

0 

3 

PIW, (mean) 

0.96 

0 

12.1 

0.54 

0.0 

9.4 

-0.54 

0.0 

9.4 


The regression coefficients c ^d, c are obtained 
by adding noise terms to linear regression coeffi¬ 
cient from experimental observations [2]: 

= &d,c T Cd, c/20, 

^d,c = Pd,c + Cd,c/2, (7) 

where a^, c , /^, c are regression slope and offset terms 
obtained from empirical data [2], The noise terms 
(d,c are taken from a normal distribution with the 
experimental value of variance a^ c obtained from 
the field data [2] . 


8.2. Equations of motion of a drifter 

The equation of motion of a floating body in the 
sea is derived from the force balance equation m : 

{rn 171 ) — = F W ave T F w ind T F ocean , (8) 

where m and rn! are masses of drifter and an ad¬ 
ditive mass, Vdr is a drift velocity , F w i n d,F wave 
and Focean are wind, wave and water drags corre¬ 
spondingly. 

Since observations show that ships and small 
size floaters da [ 23 ] reach terminal velocity rapidly 
an infinite acceleration and constant velocity during 
every time step is assumed. When the wave force is 
neglected mm Eq.[8]with implicit expressions for 
wind F wind = CaPa^UTViollWio and water drag 
Focean = C 0 p 0 A 2 0 1|' V L - dr/dt|| (V L - dr/dt) takes 
the form m- 


0 


CaPaA a \\W 10 \\W 10 

dr 


C 0 PoA 0 


V L - 


dt 


V L 


dr 

dt 


(9) 


where Vl is the Lagrangian velocity of drifter, ^ 
is the drifter velocity, C a , 0 are the drag coefficients, 
p a?G are densities and A a , 0 are effective cross sec¬ 
tional areas of body exposed to air and water corre¬ 
spondingly. Thus the equation of motion of drifters 
can be written as follows m- 

^ = V L + aW 10 , (10) 

where the coefficient a = yf p a A a C a /p 0 A 0 C 0 de¬ 
pends on the water/air constants. By expressing 
the Lagrangian velocity in terms of the Eulerian 
current and the Stockes drift Vst we obtain the 
equation [24] : 


= V e +/3W 10 , (11) 

dt 

where /? = a + ||V 5 t||/||Wio|| is derived from the 
assumption that the wind direction and the Stockes 
drift are parallel |24| . The Leeway model uses ex¬ 
plicit substitution of the second term in ( p~Q| ) by the 
Leeway drift parametrized by wind ©• 


8.3. List of categories of drifters and Leeway coef¬ 
ficients 

The Leeway coefficients for the categories used 
in the experiments are provided here. The slope, 
offset and std for the CWL, DWL are given in Ta¬ 
ble I2J 
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